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Abstract 

We present the first comprehensive fitting formula for exchange reactions of arbitrary mass 
ratios. In a comparison with numerical results, this expression is shown to be accurate 
in the hard binary limit to within 25% for most mass ratios. The result will be useful in 
forming quantitative estimates for the branching ratios of various exchange reactions in 
astrophysical applications. For example, it can be used to construct quantitative formation 
scenarios for unusual objects in globular clusters, such as binaries containing a pulsar. 
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1. Introduction 

In stellar dynamics, an exchange reaction is a particular type of interaction between 
a binary star and an incoming third star, where one of the components of the binary is 
expelled and its place is taken by the incomer. 

In the long history of the three-body problem, the study of exchange reactions had a 
curious start. Though numerical examples were given by Becker (1920) long ago, it appears 
that the very possibility of such reactions remained controversial until later numerical work 
published in 1975 (Marchal 1990). This is all the more curious because in that very year 
Hills (1975) and Heggie (1975) separately published different treatments which, taking the 
existence of exchange for granted, attempted to determine how its probability depended 
on such parameters as the initial speed of the incomer, and the masses of the stars. 

Since that time a considerable number of studies of exchange reactions have been 
carried out. In this paper we concentrate on the problem of determining the cross section 
for exchange, in the limit in which the initial speed of the incoming star is very low, but for 
all possible masses. Our study will therefore be largely complementary to previous work, 
which has often adopted different restrictions, e.g. encounters at zero impact parameter 
(Hills & Fullerton 1980), the case in which one component of the binary has very low mass 
(Hills & Dissly 1989), or the initial eccentricity is zero (Hills 1991, 1992). Other work of 
this kind will be mentioned in §5.1, for comparison with our own data. 

This paper is the seventh in a series discussing many aspects of three-body scattering 
in the point mass approximation (Hut & Bahcall 1983, Hut 1983, 1993, Heggie & Hut 
1993, Goodman & Hut 1993, McMillan & Hut 1996), but this is the first to deal with stars 
of unequal mass. For the case of equal masses, much further information on exchange cross 
sections will be found in earlier papers of this series, especially Papers II and IV, as well 
as the atlas of hard binary scattering cross sections provided by Hut (1984). 

The present paper is arranged as follows. In the following section we describe the 
numerical software we have used for generating cross sections. Because of its highly au- 
tomated yet flexible construction, this is a topic of interest in its own right. In §3 we 
analyse the problem analytically, in order to understand the dependence of the exchange 
cross section on the masses, especially in various asymptotic regimes. It turns out to be 
possible to write down a single expression which accommodates all asymptotic regimes. 
Section 4 synthesises all our numerical data and asymptotic theory to provide a compre- 
hensive and simple formula which is believed to be approximately valid for all masses. In 
the concluding section it is tested against previous work by Sigurdsson & Phinney (1993), 
Hills (1992), and Rappaport et al. (1989), and then we summarise our results. 

2. Numerical Computation of Exchange Cross Sections 

2.1 Software for Three-Body Scattering 
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The first sets of binary — single-star scattering experiments were reported by Hills 
(1975) and Heggie (1975). In the former, most encounters took place at zero impact 
parameter. The first direct determination of accurate cross sections and reaction rates for 
binary — single-star scattering was made by Hut & Bahcall (1983). For each type of total 
or differential cross section, a detailed search of impact parameter space was performed as 
a pilot study, before production runs were started. The problem with the choice of impact 
parameter (lateral offset from a head-on collision, as measured at infinity) is this: allowing 
too large an impact parameter can imply a large waste of computer time on uninteresting 
orbits; while choosing too small an impact parameter will yield a systematic underestimate 
of some cross sections, since some encounters of interest will be missed. 

The first automatic determinations of cross sections and reaction rates for binary— 
single-star scattering are described by McMillan & Hut (1996; hereafter referred to as Paper 
VI). Rather than relying on human inspection of pilot calculations, their software package 
includes an automatic feedback system that ensures near-optimal coverage of parameter 
space while guaranteeing completeness. We refer the interested reader to Paper VI for 
further details on the STARLAB software package. References to earlier papers on 3-body 
scattering can be found in the recent papers by Hills (1992), Heggie & Hut (1993), Hut 
(1993) and Sigurdsson & Phinney (1993). 

2.2 Numerical Results 

First we explain our notation. Let m\ and m 2 be the masses of the components of 
the binary, and ms that of the incoming third star. Then we define Mi 2 = m\ + 777-2, 
M123 = mi + m 2 + 7713, etc. When the incoming third body and the binary are still at 
a very large distance, let the semi-major axis of the binary be a and the relative speed 
of the third body and the barycentre of the binary be V. Then we can scale V by the 
critical value, V c , for which the total energy of the triple system, in the rest-frame of its 
barycentre, is zero, i.e. let 



All of our runs have been carried out with v = 0.1. 

We also scale the cross sections themselves, in two ways. Let S be the cross section 
for an exchange process. First, following Paper I in this series, we may scale out the 
gravitational focusing of the third body as it approaches the binary, and also the physical 
cross section of the binary itself, by defining the dimensionless cross section 



v = V/V c 



(1) 



where 



y c 2 = Gmim 2 M 12 3/(M 12 m 3 a). 



(2) 



v 



2 E 



(3) 



a = 
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This definition has the disadvantage that a becomes very large when m 3 3> M 12 , and we 
have found it convenient to scale V not by V c but by a typical speed reached by the third 
body when it makes a close approach to the binary. To be precise, we define a speed V g 

by Vg = — - 123 , and let v = V/V g . Evidently V g is the relative speed of the third body 
if it falls from rest at infinity to a distance 4a from a body of mass M 12 . We have chosen 
this slightly odd numerical factor so that V g = V c in the case of equal masses. Thus we 
define a new dimensionless cross section a = V 2 T,/(na 2 Vg), i.e. 

2V 2 E 
7rGMi23a' 

which is related to the previous definition by a = av 2 /v 2 . 

The results of our runs are displayed in Table 1, with estimates of 1-cr errors. In general 
we have attempted to ensure that the total exchange cross section for each component of 
the binary is calculated to better than about ten percent, though there are clearly cases in 
which the cross section must be so small that the resulting computational effort would be 
prohibitive. This is also illustrated in Fig.l, whose main purpose is to show the coverage 
of the plane of mass ratios in our numerical experiments. We covered all mass ratios up to 
a maximum of 2dex (between any pair of stars in the system), in steps of 0.5dex, as well 
as a few other cases. 

It is helpful to distinguish two different kinds of exchange reaction (e.g. Heggie 
1975), and Table 1 presents results for both. In direct exchange the encounter termi- 
nates promptly, and the orbits are uncomplicated; while in resonant exchange the three 
bodies form a temporary bound system, and the escaping particle emerges only after sev- 
eral interactions. The distinction between direct and resonant encounters is not always a 
clear one; the operational procedure used to classify our numerical results is presented in 
Paper VI. 

A few other remarks about the data in Table 1 should be made at this point. First, for 
a few mass ratios we have data from additional runs which are not shown here. Those shown 
are the data sets with the smallest errors. The additional data has been used, however, 
in the parameter fitting in §4 below. Second, entries with a star indicate experiments in 
which no events of the relevant kind were observed. Though it might be thought that it 
should be possible, on the basis of a large number of scattering experiments, to give an 
upper bound for the cross section of a process which produced no events, this is actually 
not rigorously possible because of the organisation of the software, which decides on the 
range of impact parameters on the basis of observed events. If the software has no evidence 
on the range of impact parameters which can produce a given event then it is possible that 
it can miss a range where the process is important. Therefore this null data should simply 
be discarded. 
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3. Asymptotic Theory of Exchange Cross Sections 

In this section we address theoretically the problem of determining cross sections for 
the reactions discussed in this paper, i.e. exchange for hard binaries. Our aim will be to 
determine the way in which they scale with the parameters of the problem, especially in the 
extreme regimes of masses. Some of these extreme cases might seem physically implausible 
or unimportant, but the purpose of the theory is to try to account for trends in the 
numerical data, and to suggest ways in which the numerical results might be extrapolated. 

In general we consider the approach of a third body whose speed "at infinity" relative 
to the barycentre of the binary was V, and denote the initial semi- major axis of the 
binary by a. We shall label the components of the binary such that mi is the mass of 
the component which is ejected as a result of the encounter, and in this section we shall 
generally add a subscript 1 to symbols for the cross section, in order to reinforce this 
convention. It must always be borne in mind that we are dealing with the case in which 
V is very small. 

3.1 The Case of a Massive Incoming Star. 

First we consider the regime in which both mi <C 7773 and <C 777.3. I n this case 
the incomer is very massive, and we shall see that it is possible for a tidal encounter by 
the third body to unbind the binary. Since we are always considering the hard binary 
limit (i.e. V <C V c ), the three stars cannot escape singly to infinity, and so an exchange 
interaction must occur. What is less obvious is to decide which component escapes. 

Let R p 3> a denote the distance of closest approach of the third body. (The subscript 
denotes pericentre.) At this distance its speed relative to the barycentre of the binary is 
denoted by V p and can be estimated from 

V p 2 ~ GM 123 /R p . (5) 

The duration of the encounter is of order R p /V p . During this period the (tidal) acceleration 
of the relative motion of the binary components by the third body is of order Gm 3 a/ R p . 
Therefore the change in their relative speed is of order AV12 ~ Gm 3 a/(RpV p ). This is 
enough to disrupt the binary if (AV12) 2 ^ GM\ 2 /a, i.e. if 

R 3 p & mla 3 /(M 123 M 12 ). (6) 

Incidentally we have assumed implicitly in all this that the duration of the encounter does 
not much exceed the period of the binary, for then the change in energy of the binary 
would be exponentially smaller than the estimate we have used (Heggie 1975). In fact it 
is easy to show that the assumption on time scales is justified post hoc by eq.(6). 
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If, as we assume throughout, the original speed of approach (V) of the third body 
is very small, the barycentre of the binary moves along a nearly parabolic orbit relative 
to 777-3. As a result of the disruption of the binary, one component will now be moving 
in front of, and faster than, the barycentre, while the other will fall behind, and the 
probability that a given component is the one that moves in front is roughly the same for 
both components. Call this body nif and the other component nib- If m/ nib, then nif 
is moving only slightly faster than the barycentre, but since the barycentre is moving on a 
slightly hyperbolic orbit, it follows that nif will escape. If, on the other hand, m/ nib, 
then nif will pick up almost all the extra velocity (AV12) generated in the encounter, which 
is enough to ensure its escape. In either case, then, m/ is the component of the original 
binary which escapes, and, by the convention already stated, we identify this with ni\. 

It might be surprising to find that the probability that a given component is the escaper 
does not diminish significantly when it is much more massive than the other component. 
This is, however, confirmed in Table 1 by such cases as ni\ = 0.09, m 3 = 9.09. 

It remains only to estimate the cross section for this process. The impact parameter 
of the third body, p, when at a large distance on its incoming orbit, is related to other 
parameters of the encounter by approximate angular momentum conservation in its motion 
relative to the barycentre of the binary, which leads to pV = R p V p . Therefore the cross 
section for these encounters is Si ~ p 2 ~ (R P V P /V) 2 . (The factor (V p /V) 2 accounts for 
the "gravitational focusing" that takes place as the third body approaches on its nearly 
parabolic orbit.) Using eq.(5), and the expression on the right of eq.(6) to estimate R^, 



In view of our assumptions that ni\ and 777,2 are very small compared with 7773, we see 
that an equally valid asymptotic expression is 



This will be more useful for combining with other asymptotic formulae later. 

Finally we have to dispose of the question of resonant exchange. In fact, if, after this 
first encounter, the third body is bound to the binary without disrupting it, thus forming a 
resonance, then on each subsequent encounter it is quite likely to become unbound again, 
and so the cross section for resonant exchange cannot greatly exceed that just estimated 
for direct exchange. Again this is confirmed by data in Table 1, though it is noticeable 
that the cross section for ejection is now somewhat larger for the less massive star. 

This discussion has been taken fairly slowly so as to introduce the various steps and 
considerations involved in the estimate of the exchange cross section in a certain regime. In 
subsequent discussions only the distinctively different aspects will be treated in comparable 
detail. 





(7) 



6 



There are several aspects of this kind of encounter which will be useful later (§3.2). 
The energy given to the binary will be comparable to the energy required to break it up. 
Therefore the energy extracted from the relative motion of 7773 and the barycentre of the 
binary must also be comparable with Gmxm^i 'a, and so this can be taken as an estimate 
of the binding energy of the new binary. It follows that the semi-major axis of the new 
orbit, denoted by 023, may be estimated by 

Gmim2/a ^ Gmsm 2 /a23, (8) 

i.e. 023 ~ 07713/777-1. Also the pericentric distance of the new binary is of order R p , and so 
can be expressed as 

R P ~ a 2 3(m 1 /m 3 )(M 1 23/M 12 ) 1/3 . (9) 

It follows that the eccentricity of the new binary (e') is high, since 1 — e' ~ R p /ci23- Finally, 
the initial separation of mi and 7712 is related to the new semi-major axis by 

R12/CI23 ~ K/m 3 ). (10) 

The argument that the binding energy of the new binary must be comparable with 
that of the original one has analogues in other capture processes. One familiar example 
is tidal capture in globular star clusters (Fabian et al. 1975), where the encounter must 
be close enough to remove the kinetic energy, T, of relative motion of the two stars. In 
this case it follows that the binding energy of the new binary will be comparable to T. 
The analogy is not precise, however, because in this case the binding is not achieved by 
disruption of one of the stars. 

3.2 The Case of Ejection of a Massive Star. 

We move on now to the regime in which 7772 <C mi and 7773 <C mi The first point to 
notice about this case is that it may be obtained from the previous case by time-reversal. 
In the previous case a heavy incoming third body becomes a component in the new binary, 
while in the present case it is the heavy component which is ejected. (Recall that the 
components of the binary are named in such a way that it is the first component which is 
exchanged.) While the ejection of the most massive component in a three-body interaction 
seems unlikely, the use of time-reversal yields the special circumstances in which it can 
happen. In fact the encounter must occur while the two components of the binary are close 
to pericentre (since R p <C C123 in the previous case) , and the incoming star must come even 
closer to the lighter component, 7772. 

We shall return to some aspects of a direct calculation of the cross section towards 
the end of this section. Mainly, however, we shall exploit the fact that time-reversal leads 
to a process whose cross section we have already estimated, and so we may use the theory 
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of detailed balance to estimate the required result. This theory is set out in the Appendix, 
and leads to the following result for inverse process. (Here we omit the subscript 1 on the 
cross section, since the identity of the escaper is sufficiently defined by other aspects of the 
notation.) 

Let — — — (E 12 — > £23) be the differential cross section for the formation of a binary 

having components m2, and energy £ 23 , from a binary having components mi, m 2 
and energy £12, by exchange. Then this is related to the differential cross section for the 
inverse process by 



dE JEi 3 ^E l ,)=\ K -J j y,{^) ^(fe - E 23 ). (11) 

Here Vi and V3 are the speeds of the incoming single body in the two scatterings. The 
coefficient on the right side comes from the phase space volumes associated with the re- 
actants: for the incoming single body this is proportional to V 2 , and the factor E~^ 2 is 
easily understood in a similar way in terms of the internal degrees of freedom of the binary 
(Hut 1985). 

This relation involves differential cross sections, whereas eq.(7) estimates a total cross 
section. However, we already estimated that the typical energy of the new binary in that 

case is £ 23 ~ £i 2 , and if we estimate (£23 — »■ £12) ~ £(£23 ~~ * -Ei2)/£i2, we see that 

a£i2 

eq.(ll) yields the following estimate for the total cross section: 



E(£ 23 -> £12) ~ ( ^S(£ 12 -> £ 23 ). 



Substituting eq.(7) in the right side we find, for the case mi, m 2 <C m 3 , the estimate 

S(£ 2 3 — > £i 2 ) 



GM 123 a 12 fM 123 \ 1/3 /mi \ 5/2 f 1/2 



V{ V M 12 J \m 3 J \M 23/ 

where ai2 (heretofore denoted simply by a) is the semi-major axis of the binary with 
components 1 and 2. Now we use eq.(8) to replace ai2 by a23mi/m3, and interchange the 
labeling of stars 1 and 2 to restore our customary labeling of the incomer. This yields the 
exchange cross section 

V3 ,^7/2 / M23 XV2 



M12, 

in the case m-2, m 3 <C mi. An equally valid asymptotic formula is obtained by replacing 
Mi 2 by M123, and then we see that a formula which is compatible with eq.(7), and therefore 
valid in both the asymptotic regimes of §§3.1 and 3.2, is 

^ GMrxa / M^ ^ / m^ ^ ( ^ 

V 2 \M 123 ) \M 13 J \M 12 J 1 ' 
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Before we leave this regime it is worth noting that the major part of the mass- 
dependence here, i.e. the factor (m 3 /M 13 ) 7 / 2 is easily understood. We already saw (eq.(9), 
but with relabeling appropriate to time-reversal) that the separation of the binary com- 
ponents at the time of approach of m 3 must be of order R12 ~ a(m 3 / 'mi) (M "123 / 'M23) 1 "' /3 . 
Now the probability of this (for a thermal distribution of binaries of a given energy) is of 
order (i?i2/o) 5 / 2 when R12 <C a (cf. Paper IV). Next, the distance of closest approach 
of m 3 to ni2 must be of order R p ~ arris/ mi, by eq.(10), and at this time the speed 
of the third body, which is gained mostly by falling to within a distance R12 of mi, is 



ft2y2 / R ^5/2 

given by V 2 ~ GMi2s/Ri2- It follows that the cross section is Si ~ ( — 

V 2 \ a 

7 /2 \ 5/6 



~ = — — . This is larger than our estimate in eq. 12 because we 

V 2 \m 1 J \M23J 

have not taken into account the special circumstances of the encounter which allows m 3 
to be captured by m-2, but does explain the major part of the mass dependence: all three 
stars must come within a separation which is of order the pericentric distance of a binary of 
high eccentricity. The factor (m 3 /mi) 7 / 2 is a measure of the phase space volume available 
to such systems. 

3.3 The Case in Which a Massive Component Remains 

Finally we turn to the regime in which both mi <C m2 and m 3 <C mi. We first 
discuss the case in which m 3 > mi, i.e. an object of intermediate mass displaces a low- 
mass companion of an object of very high mass. (Recall our convention for labelling the 
components, which is that mi is the mass of the component which is ejected.) To begin, 
let us consider the possibility of direct exchange. Suppose m 3 approaches mi within a 
distance R p <C a. Then the speed of m 3 will be given by V 2 ~ GMi2s/a provided that the 
influence of mi is not dominant, i.e. provided that R p J> C1M13/M123. It follows that the 
speed imparted to mi is of order the escape speed provided that R p ~ aMi 3 /Mi 2 3. Hence 
the cross section for direct exchange is 



/ \ 2 

GMi 23 a / m 3 



V 2 \ TO 2 



(13) 



Now suppose only that the distance of closest approach between m 3 and mi is of order 
a. Then there is a significant probability that m 3 will become bound to the binary without 
ejecting mi, i.e. a resonance will form. Thus the cross section for resonance formation 
greatly exceeds the cross section for direct exchange, since our estimate of R p for direct 
exchange requires that R p a. In fact the cross section for formation of a resonance is 
simply 

GMi 23 q 

Si ^ 2 —. (14) 
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At this point in the evolution of the resonance the binding energy of ms is of order 
the change in energy of mi, which is of order Gm\m^/a. Note that this estimate is valid 
no matter how small mi is; we are always considering the limit of extreme hardness, and 
even a small change in the energy of the third body can bind it to the binary if its energy 
at infinity was sufficiently small. On the other hand, the cross section for formation of a 
resonance leading to exchange cannot greatly exceed our estimate: if the closest distance 
of approach of ms greatly exceeds a then a hierarchical triple system forms, and exchange 
is very improbable. 

Now we must estimate the probability that the resonance will be resolved with the 
escape of star m±. We think of the binding energy of this particle, E\, as performing a 
random walk under the influence of repeated passages by star ms- The typical change in 
Ei is of order Gmsmi/a. Because we are assuming that ms mi, it might be thought 
that the mean effect would be to systematically unbind mi, by a kind of mass segregation. 
However, the mean change in E\, taken over an ensemble of such systems, is actually of 
second order in the ratio of the masses: in the approximation of first-order perturbation 
theory, time-reversal shows that for each change in Ei there is a system in which the 
change has the opposite sign. Therefore we shall assume that the mean change for a given 
system is negligible. 

As an aside, it is worth mentioning here that the system is a kind of hierarchical triple. 
Usually this term is used in reference to a binary about which a third star revolves on a 
large elliptical orbit which is well separated spatially from the binary. In that case the 
perturbation of the third body is weak because the orbit of the third body is large. In the 
present case there is no such spatial separation, but still the perturbation of the third body 
by the binary is weak, and this is due to the low mass of one component of the binary. 

We must estimate the probability that Ei, starting from the value of order Gmim2/a, 
may randomly walk, in steps of order Gmims/a, to the value without first reaching the 
value Gmim2 / a + Gmims / a (as ms would then escape again). Now the probability that a 
one-dimensional random walk exits from a given boundary is a linear function of the initial 
position, varying from unity at the boundary of interest to at the opposite boundary. 
In the present case the boundaries are at E\ = Gmim2/a + Gmims/a and 0, and so the 
probability of escape at the boundary E\ = 0, starting at E\ = Gmim^/a, is of order 
m 3 /m 2 . This can be estimated equally well as -M13/M123, and so it follows, using eq.(14), 
that the cross section for resonant exchange is 



It also follows that we obtain a form which is asymptotically correct in all regimes studied 
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so far if we modify eq.(12) to 



GM 123 a ( M 23 \ 1/6 / m 3 \ 7/2 / M 123 \ 1/S ( M 13 



Sl ~ v- vm 123 J Im 13 J Im 12 J Im 123 J- (l5) 



or, if everything is normalised by the total mass of the triple system, 
Si 



GM 123 a ( M 23 \ 1/6 f ™ 3 \ 7/2 f M 123 \ 1/3 ( M 13 ^ 5/2 



V 2 \M 123 J \M 123 J \M 12 J \M 123 , 

We also note that the energy of the new binary will be comparable with that of the initial 
binary. 

Observe that our estimate for resonant exchange is indeed larger than the estimate, 
eq.(13), for direct exchange. The fact that this is true in the mass regime under discussion 
is also illustrated in Table 1 by such cases as mi = 0.01, m 3 = 0.099. 

Finally we turn to the case m 3 <C mi <C m 2 , which is the time reversal of the case 
just considered. Though the foregoing argument for the formation of a resonance goes 
through, it is now very unlikely to be resolved by exchange. We can, however, estimate 
the cross section for resonant exchange by detailed balance, using the same method as in 
§3.2. The result is that 

V2 \ 1/2 



GM 123 a (rnzY 1 ( 1 ( M 13 \ 
V 2 \mj \M 12 J \M 123 J 



and we see that eq.(15) is, once again, a valid asymptotic result in this regime. We therefore 
adopt eq.(15) as a cross section whose form is valid in all regimes. Henceforth we drop 
the subscript 1 on E, and remind the reader that m\ is the mass of the component of the 
original binary which is ejected. 

4. Synthesis of Numerical and Analytical Results 

4.1 A Test of the Asymptotic Formulae 

Though we already made a number of qualitative remarks relating the theory of §3 to 
the numerical data in Table 1, it is now time for a more quantitative study. Our aim is not 
yet to provide a fitting formula for the data (which we take up in §4.2), but nevertheless 
we shall generalise eq.(15) slightly. Using the notation of eq.(4), we may rewrite eq.(15) as 

, -1/6 / \ -7/2 / \ 1/3 / N -1 



M 23 J V m 3 J \ M 12 J V M 13i 

In order to allow different multiplicative constants in the different asymptotic regimes, we 
generalise this to 

/ \ l/ 3 / \ -V 2 / \ "I/ 6 / \ -1 

a ~ ai a 2 + —— a 3 H o 4 + — — a 5 + — — - , (16) 

V ^12/ V ^3/ V M 23j V M is) 
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where the ai are constants. There is then one asymptotic form for the regime of §3.1, 
one for that of §3.2, and one each for the two regimes in §3.3 (i.e. the regimes in which 
mi <C m 3 and mi ^> m 3 ). 

The formula in eq.(16) has been fitted to the logarithm of the data in Table 1 using 
least squares, the standard deviation in the logarithm being estimated by the relative error. 
In cases where no cross section was measurable the data point was assigned zero weight. 

Since the aim of this exercise was to detect possible errors in the theoretical asymptotic 
form we searched for deviations (between the results of this formula and the experimental 
data) which were significant ("2-cr"), large (above 1 in the natural logarithm), and in an 
extreme mass regime. We illustrate the results by considering in a little detail the set of 
data in which the discrepancies were largest. This was the case m 2 <C m 3 ~ mi/3, as 
illustrated by Table 2. The numerical results show a trend, but not a significant one, and 
a constant value in the numerical data (as predicted by the formula) is not ruled out at the 
20% level. A very similar set of results, but with smaller discrepancies and in the opposite 
sense, is found in the case m 2 « mi = m 3 . Apart from these, the only systematic, large 
discrepancies occur at one or two data points where the masses are comparable, and the 
asymptotic formulae need not apply. In conclusion, then, there is no evidence that the 
theoretical formula is inconsistent, in the appropriate regimes, with trends in the numerical 
data. 



4.2 A Semi-Numerical Fitting Formula 

The theoretical results of §3 are intended to provide the asymptotic dependence of 
the exchange cross section on the masses of the participants, but do not even attempt to 
provide numerical coefficients for these. The numerical data, on the other hand, apply to 
only discrete points in the parameter space of the masses. An obvious way of synthesising 
the two kinds of result is to adopt a form with the same asymptotic properties as the 
analytical result, but with additional terms which are chosen to optimise the agreement 
with the numerical data. In a sense this is what was done in the previous section, but 
with a limited degree of flexibility. Here we adopt a more general approach which could 
be extended more or less arbitrarily. 

Again we switch to a, defined by eq.(4), and generalise eq.(15) to 



a = 



Jl/f \ V 6 / \ V 2 / Jl/T \ 1/3 / »i \ / rn+n=N 

-M23 \ 1 m 3 \ I M123 \ / M- - • 



M 123 J \M 13 J \M 12 
where the a mn are constants, 



\m,n;m+n=0 ' 



fj,i = mi/Mi2, and fi 2 = m 3 /M 12 3- 
12 



(18) 



These two parameters entirely span the possible ranges of mass ratios in the unit square 
< [ii < 1. The exponential is used to constrain the function to be positive, and 
also to avoid altering the asymptotic character of the leading expression. By taking larger 
values of the highest power N, it would be possible, in principle, to improve the fit to 
arbitrary accuracy. 

We have fitted formulae with N < 5, and at the largest value the value of x 2 is 133, 
which, considering the number of degrees of freedom, is still rather large. (The number of 
data points is 126, and the number of free coefficients is 21.) Nevertheless we suggest the 
use of a cubic polynomial in the exponential in eq.(17), with the coefficients given in Table 
3. As these 10 coefficients are quoted to 2 decimal places, the maximum relative error in 
the evaluation of a is about 5%. This is adequate in view of the accuracy of the fit, which 
is discussed further below. 

One worry about using a single polynomial for any kind of interpolation is the possi- 
bility of large oscillations between data points, but in fact the polynomial with the above 
coefficients is well behaved. Its range is about 5, i.e a varies by a factor of about 100 by 
the effect of this polynomial. Its most noticeable feature is a minimum value at [ii = 1, 
[t2 = 0, i.e. the regime mi ^> M23. The cross section is very small in this region anyway, 
because of the factor (m 3 /M 13 ) 7 / 2 in eq.(17). This can also be seen in Fig.l, which shows 
contours of log a. 

The fitting formula is quite successful. For half of our measurements the result is 
accurate to better than 10%, and for about 75% of our measurements it is better than 
20%. Of the remaining measurements there are some in which the disagreement exceeds 
2 standard deviations, and they are shown in Fig. 2. At each point the label gives the 
relative error in the sense (&th — &num) /&th, where ath and a num are the theoretical value 
(i.e. from eq.(17)) and the numerical value (from Table 1), respectively. Thus a value of 
— 1 would mean that a t h = 0.5a num , while the extreme positive value of 0.69 means that 
ath — 3a nurn . Note that the cross sections around this point are very low (Fig.l), which 
has two consequences: first, that the standard deviation of the numerical cross section is 
comparable with the cross section itself, and, second, that exchange may be unimportant 
in this regime. Indeed for all but two of the points plotted in Fig. 2 the discrepancy 
between the formula and the numerical data is less than three standard deviations. The 
two exceptions are the points labeled 0.53 and —0.75. It is also reassuring that almost all 
the discrepant points are surrounded by data points (Fig.l) where the agreement between 
the numerical and semi-theoretical results is satisfactory by the above criteria. Thus 
the errors are localised. Nevertheless, it is evident that there are some significant and 
systematic trends in this data, and the possible effect of these discrepancies should be 
assessed in any application of the fitting formula. 

5. Discussion and Conclusions 
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5.1 Comparison with Other Authors 

Before summarising our findings, our immediate task will be a quantitative comparison 
of our fitting formula with some existing results in the literature. The first of these that 
we shall examine is the paper of Sigurdsson & Phinney (1993). They give results for 
different speeds V, and we have chosen the data for the lowest speed, since the validity of 
our conclusions is restricted to the case of hard binaries. Data in their Tables 3A and 3B 
have been normalised, where necessary, to the cross section a (eq.(3)), and collected in our 
Table 4. Also included are results of our fitting formula, eq.(17), again converted to a. 

Unfortunately Sigurdsson & Phinney do not give estimates of the errors of their results, 
and in most of their runs the initial eccentricity of the binary was chosen to be zero, and 
so only an informal comparison is possible. The agreement is often quite good, and almost 
always within a factor of two. In the cases where the disagreement is most serious it 
is probably attributable to the different choice of initial eccentricity distribution. For 
example, in the equal-mass case our fitting formula agrees with our numerical data to 
within 10%. In the case <C m\ = m^, however, where their result for ejection of the 
low-mass component falls below that of our fitting formula, Sigurdsson & Phinney have 
missed a significant fraction of exchange encounters by too small a choice of the maximum 
impact parameter (E.S. Phinney, pers. comm.). The fact that each of their data points is 
a weighted average over a range of speeds V may complicate the comparison further. 

Now we turn to data presented by Hills (1992) on the case in which mi = m 2 . The 
initial eccentricity in his experiments was again 0, but we find that our results agree with 
Hills' to within about a factor of two over the entire range for which he found exchange 
events, i.e. 0.3 ^m 3 /mi <> 10 4 . Hills' result exceeds ours except at the lowest mass ratio 
in this range, and the discussion at the end of §3.2 suggests that the different choices of 
initial eccentricity provide a plausible explanation for this last fact. 

Finally in this section we compare our results with those of Rappaport et al. (1989). 
They computed the exchange cross section, by numerical scattering experiments, for a 
sequence of binary pulsars. The sequence is characterised either by the orbital period or 
the mass of the low-mass companion of the neutron star. The cross section was computed 
for an environment containing a stellar population drawn from dynamical models of two 
globular clusters. Here we restrict attention to their binary of shortest period (3 days), 
since our results are restricted to the hard-binary regime, and to the model of u Cen, 
for illustration. The results of Rappaport et al. give the dimensionless scattering cross 
sections T,/(na 2 ) = 1.3 and 55, for ejection of the neutron star and low-mass companion, 
respectively. Typical uncertainties are about 40% and 15%, respectively. 

For our comparison we have used eq.(17) for each of the ten components in the stellar 
population listed by Rappaport et al, and have summed the contributions, account being 
taken of their relative number density and velocity dispersions. Expressed in terms of the 
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quantity computed by them, our results are £/(7ra 2 ) = 1.37 and 76.0. The agreement 
is acceptable, considering the typical errors in all results, and the fact that the result of 
Rappaport et al. applies to an initially circular binary. It also illustrates the utility of our 
results, as the cross section could be obtained for any reasonable stellar population with 
little extra work. 



5.2 Conclusions 

This paper is a contribution to the theory of three-body classical gravitational scat- 
tering. This is a large topic with an extensive literature, but it is the application to the 
dynamics of globular star clusters which has provided the focus for our work. From the 
point of view of this application, one of the most important processes is exchange, whereby 
an incoming star ejects one component of a binary and forms a new binary with the other 
component. In the context of star clusters we are also mainly concerned with encoun- 
ters with hard binaries, which are too energetic for an encounter to break up the system 
into three single stars. Finally, this application dictates the importance of understanding 
scattering in a system where the stars may have quite widely differing masses. 

The main result of this paper is a semi-analytical cross section for exchange, in the 
hard-binary limit, for all possible masses. It has been derived partly from theoretical con- 
siderations and partly from extensive new numerical data on scattering events. The theory 
allowed us to estimate the dependence of the cross section on the masses, in various limit- 
ing cases, and the numerical data showed how the theoretical results can be parameterised 
so as to provide a better fit, including the cases where the masses are comparable. 

The result is given in eq.(17) with coefficients in Table 3, and we now summarise the 
way in which this information may be used. Suppose a star of mass rris approaches a 
binary with components of mass mi, m-2, and that its speed, while still at a large distance 
from the binary, is V relative to the binary. Let the initial semi-major axis of the binary be 
a. Then the cross section for events in which the component of mass mi is ejected, leaving 
a binary consisting of the other two stars, can be computed in the following way. First 
compute M12 = mi +m2, and also M23, M13 and M123, defined similarly. Next, compute 
a from eq.(17), where N = 3, /ii and ^2 are defined in eqs.(18) and the coefficients are 
taken from Table 3. (Note that the exponential in eq.(17) is simply that of the cubic 
aoo + ttio^i + a oiA*2 + • • • + ai2/Ui//| + ao3/U2-) Then the required cross section is given by 
solving eq.(4) for E. This result is approximately valid provided that the binary is hard, 
i.e. v 2 <C 1, where v is defined by eq.(l). 
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In astrophysical units this can all be summarised in the formulae 

2 / n/r \ / n/r \ 1 / 6 / \ 7 / 2 / n/r \ 1 / 3 / n/r 
> i . i M 23 \ / m 3 \ / M123 \ / Mi 3 



X 



,Mi 23 / V M 13/ \M 12 J \M U 
x g 3. 70+7. 49^1-1-89^2 -15. 49^^ -2. 93miM2 -2. 92^2+3. 07^i + 13. 15^?M2 -5. 23^1^2+3.12^2 ^TJ 2 

(19) 

where = m\/M\2 and ^2 = "W-M123. The condition that the binary is hard is 

Oil (— V — \ 2 ( -?-) ( M ^ M A « 1 ( 2 o) 
U - Uii Vl0km/s J UlAlJ \m 1 m 2 M 123 ) * A ' ^ UJ 

(The last factor is dimensionless, and the factor M may be omitted if the masses are 
expressed in units of the solar mass.) 

For example, consider a binary pulsar consisting of a neutron star of mass 1.4M© and a 
white dwarf companion of mass 0.2M© in an orbit of period 10 days. Suppose it encounters 
a 1OM black hole with a relative speed of lOkm/s. Then a ~ 0.11AU, \i 2 — 0.862, and 
the left side of eq.(20) evaluates to 0.06. To compute the cross section for ejection of the 
neutron star we have fii = 0.875, and so eq.(19) gives E ~ 120AU 2 . Surprisingly, perhaps, 
the cross section is not much smaller than that for ejection of the low-mass companion; this 
is obtained by setting jii = 0.125, and is approximately 200AU 2 . Thus, of all encounters 
leading to exchange, roughly one third lead to ejection of the companion of higher mass, 
for these parameters. 

Some cautionary remarks are now in order. First, the cross section in eq.(19) is a 
statistical result, and one of the averages that is implicitly performed in our work is over 
the initial eccentricity of the binary. It is assumed to have the "thermal" distribution 
/(e) = 2e. The limited evidence available (§5.1) suggests that this affects the result by 
at most a factor of two, except in regimes where exchange is probable only when the 
eccentricity is high. From the discussion of §3.2 (which also applies with some changes to 
§3.3), these events are those in which <^ mi; in such circumstances exchange is actually 
very rare anyway, and so this issue is unlikely to be important. 

The second precaution concerns the error in the fitting formula. We have found that 
it usually agrees with numerical data to better than 20%, but that there are a few places 
where the error can apparently exceed 50%. These are illustrated in Fig. 2, and it may be 
advisable to check whether events in these areas of parameter space are of importance in 
a given application. 

Finally, no-one who makes use of these results in applications will need to be reminded 
that they apply to the point-mass approximation. In some cases the results would be 
drastically different for stars of finite radius. 

We thank the referee for his detailed and helpful comments. This work was supported in 
part by NASA grant NAGW-2559 and NSF grant AST-9308005. DCH thanks the Institute 



16 



for Advanced Study, Princeton, for its hospitality while much of the theoretical work in 
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Note: the columns headed "1" and "2" give the cross sections for exchange in which the 
particle of mass mi or m2, respectively, is ejected. 
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Example of Discrepant Results 
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Coefficients for a Semi-Numerical Exchange Cross Section 
m 1 2 3 
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3.70 7.49 -15.49 3.07 

1 -1.89 -2.93 13.15 

2 -2.92 -5.23 

3 3.12 
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Table 4 

Comparison with Results of Sigurdsson & Phinney for a 
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Note: the columns headed "Star 1" and "Star 2" give the cross sections for exchange in 
which the particle of mass mi or rn,2, respectively, is ejected. 
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Appendix. Detailed Balance for Exchange Cross Sections 

The theory of detailed balance is described in some generality in Heggie (1975), 
though it is expressed in terms of rate functions, i.e. the integral of a cross section over a 
Maxwellian distribution of velocities, and does not explicitly deal with exchange reactions. 
Detailed balance is also described in terms of cross sections in Paper III in this series, 
though in the case of equal masses. Since the integration over a Maxwellian is essentially 
a Laplace transform, it is possible to obtain the result for cross sections from the result 
in Heggie (1975), and this will be the starting point for the following treatment. We have 
verified that a direct derivation for cross sections for exchange reactions with different 
masses leads to the same result. 

With some changes of notation the result presented in Heggie (1975) can be written 

as 

1 dR 

-n 1 n 2 n 3 {Ti/kTf/ 2 {m 1 m 2 f\E l2 \- 5 ' 2 — r {E 12 - E[ 2 ) exp(-E 12 /kT) = 

Z &illy2, 

1 dT? 

^n 1 n 2 n 3 (7r/fcT) 3 / 2 (m 1 m 2 ) 3 |i?; 2 |- 5 / 2 -^(^ 2 -> E 12 ) exp(-^ 2 /fcT), 
Z ai±i\ 2 

(A.l) 

where dR{E\ 2 — > E' 12 ) is the rate (per unit density of reactants) of reactions which change 
the binding energy of a binary from E 12 to a value E' 12 within a range of size dE' 12 , T is the 
kinetic temperature, and n; is the number density of stars with mass m t ; these cancel from 
this equation, and are irrelevant in what follows. Eq.(A.l) is appropriate to encounters 
which do not lead to exchange, and the use of the labels 12, which identify the components 
of the binary, seems pedantic at this stage, but it becomes useful when we go on to discuss 
exchange. 

Now the rate function can be defined in terms of a differential cross section by 



dE' 12 v 1/7 V m 3 M 12 J J V 2M 123 kT J dE[ 



'12 



Then we can find the detailed balance relation for the differential cross section by substi- 



tuting this integral into eq.(A.l) and inserting a delta function 5 ( 1Z V 3 2 + E\ 2 — E 

\ 2M 12 3 

which isolates interactions involving systems with total energy E in their barycentric frame. 
(Primed variables are used on the right side of eq.(A.l).) Cancelling symmetric functions 
of the masses, the result we obtain is 

(m 1 m 2 ) s (m 3 M 12 ) 1 / 2 \E 12 \- 5 / 2 Vi——(E 12 -> E' 12 ) = 



12 T/2 



dE[ 2 

{ mi m 2 f(m 3 M 12 yi 2 \Ei 12 \-*/ 2 v£ ~^(E[ 2 - E 12 ). 



23 



Now we observe that a similar relation holds for exchange reactions, provided that 
the masses are correctly identified. Thus 

(m 1 m 2 ) 3 (m 3 M 12 ) 1 / 2 \E 12 \- 5 / 2 Vi——(E 12 -> E' 23 ) = 

"-^23 

(m 2 m 3 ) 3 (m 1 M 23 ) 1 /2|^ 3 |-5/2 y/ 2 ^_ (jE; / 3 _ ^ 

We can now drop the primes as the start and end states are sufficiently identified by the 
subscripts, and we deduce that 

cffi /rn . /mi\ 5/2 /M 12 \ 1/2 Vi /E 12 \~ 5/2 



(^23 -> E 12 ) = — 772 ( ~E< ) TFf-(^i2 -> #23), 



^i2 V 7 W/ V M 2 3 ; ^/ W dE : 

which is eq.(ll) in this paper. 



■>23 
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Figure Captions 



Figure 1: Coverage of the parameter space of mass ratios in the numerical experiments. 
Open circles represent experiments where the cross section was too small to be measurable. 
In this figure, mi is the mass of the component which is ejected. Dashed lines are contours 
of the logarithm of the theoretical exchange cross section log 10 a given by eq.(17). The 
values of log 10 a range from —5 at lower right in steps of 0.5 to 1. 

Figure 2: Data points where the fit of the semi-analytical formula, eq.(17), is relatively 
poor. At each value of the mass ratio where the relative error exceeds both 20% and two 
standard deviations, the relative error is printed. 
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